Fractal nanostructures with the Hilbert curve geometry as a SERS substrate 
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We suggest a new type of substrates for the Surface Enhanced Raman Scattering measurements 
with the geometry based on self-similar fractal space filling curves. As an example, we have studied 
theoretically the dielectric response properties of doped semiconductor nanostructures, where the 
conducting electrons are trapped in the effective potential having the geometry of the Hilbert curve. 
We have found that the system may exhibit the induced charge distribution specific for either two 
dimensional or one dimensional systems, depending on the frequency of the external applied field. 
We have demonstrated that with the increasing of the depth of the trapping potential the resonance 
of the system counterintuitively shifts to lower frequencies. 



Quick, reliable and selective single molecule detection 
of different molecular species is a challenging problem. 
The Surface Enhanced Raman Scattering (SERS) tech- 
nique is considered as one of the most powerful meth- 
ods for single molecular detection [l|, Q • The local elec- 
tromagnetic field enhancement on random metallic sub- 
strates of irregular, fractal-like geometries, or on corru- 
gated metal surfaces was considered as one of the pri- 
mary enhancement mechanisms for the SERS measure- 
ments since the discovery of the phenomenon. Such ge- 
ometry results in a local field intensity enhancement by 
a factor up to \E\'^ cx 10^ - lO'^ in so called "hot spots" 
[^, The Raman signal is proportional to \E\'^, there- 
fore a molecule in a "hot spot" may emit up to 10^ — 10^" 
times stronger signal. Such kind of substrates were in- 
tensively studied both experimentally and theoretically 
[J, 0, [B-[3| • Foi' example, it was shown, that the statistics 
of "hot spots" in such substrates may be governed by 
power law distributions 0]. However, random substrates 
have significant disadvantage of poor reproducibility of 
measurements for different substrate samples. Another 
disadvantage of these type of substrate geometries is a 
relatively low degree of controllability of the resonance 
properties of the substrate. 

In the last decade there have been a considerable 
interest in stable, reproducible SERS substrates. For 
example, arrays of nanoparticle pairs 0, Q or pho- 
tonic/plasmonic crystal structures were considered re- 
cently as an alternative geometry with more degrees of 
control of the substrate's resonance properties [10|. Such 
kind of designs can generate a periodic pattern of a high 
intensity "hot spots" , however, a relatively large fraction 
of the area is SERS inactive. In this paper we consider a 
new type of geometry for SERS substrates, which com- 
bines the advantages of irregular fractal substrates, re- 
sulting in a strong local field enhancement, and periodic 
plasmonic structures, characterized by relatively easy re- 
producibility and a high degree of controllability of the 
substrate's resonance properties. We suggest to use for 
the SERS substrate a geometry derived from a special 
class of fractal-like space filling curves llj . For example, 
the Hilbert curve is a special case of space filling curves. 



It has a large number of turns and corners for a given 
length of a curve (please, see Figs. (1,2)). Note, that the 
Hilbert curve is a self similar structure, which may be 
built making a finite (or infinite) number of iterations. 
A nanoscale SERS substrate will not be a true fractal 
because of the lack of an infinitely fine spatial resolution. 
In particular, at the atomic scale the quantization of elec- 
tronic states will start to play a significant role. Since the 
local field enhancement is the most significant at surfaces 
with a large local curvature (at the corners), the Hilbert 
curve is a good candidate as a substrate for the SERS 
detection with a high probability for a molecule to hit 
a "hot spot. Another advantage of the Hilbert curve is 
that while it has a self-similar fractal structure, its Haus- 
dorff dimension is 2, i.e. the Hilbert curve uniformly 
covers surface without empty spots. Regular fractal-like 
structures have been considered for nano-photonic ap- 
plications before l3"15|. However, in these works the 



authors considered fractal structures with the Hausdorff 
dimension 1 < D < 2, and therefore, such structures do 
not fill surface uniformly. 

Space filling substrates can be manufactured using 
modern lithography techniques [l3] ■ A possible approach 
is to use ion beam lithography methods to etch trenches 
with the shape of the Hilbert curve on a metal surface 
or a thin film. In a case of a thin film one may vary 
the depth of the trenches, continuously moving from two 
dimensional geometry of a thin film to the one dimen- 
sional self-similar geometry of a Hilbert wire. Another 
alternative approach is to use metal nanospheres orga- 
nized on a surface as the Hilbert curve. In the last case 
one has an additional degrees of control of the resonance 
properties. By choosing certain spheres radii and the 
inter-sphere distance one may shift the resonances in the 
most suitable range for the single molecule detection of 
a given chemical specie. In order to increase the de- 
gree of control of the resonant properties of the system 
one may fabricate the nanostructure using spheres made 
of different materials (for example, silver and gold), or 
using bimetal spheres [l6l |. One i xiay also use Surface 
Tunneling Microscope techniques [17| to assemble and 
measure the response of an atomic scale structure hav- 
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ing the geometry of the Hilbert curve. In order to utihze 
rich dielectric properties of fractal-hke space filhng curves 
in modern nanoelectronics one may use photo lithogra- 
phy methods, similar to those for microchip fabrication, 
which will result in non-uniform doping levels in a semi- 
conductor substrate. Conducting electrons will be local- 
ized in a highly non-uniform trapping potential with the 
Hilbert curve geometry. Note that while the electrons 
will be trapped in a quasi-lD potential, the interaction 
with an external field will result in the induced charge 
distribution having also quasi one dimensional proper- 
ties. However, the electrons may also interact with each 
other via the Coulomb interaction, which will extend out 
of the plane of the substrate, making the system effec- 
tively two dimensional. As we will show in this paper, 
the collective excitations in such systems may exhibit the 
induced charge distributions, which is characteristic for 
either one or two dimensional systems, depending on the 
excitation frequency. 

In order to calculate the dielectric properties of frac- 
tal nanostructures we used linear response theory un- 
der Random Phase Approximation (RPA) [l^]. This ap- 
proach allows us to take into account the non-local prop- 
erties of the dielectric response of the highly inhomoge- 
neous atomic scale system. 

The Schrodinger equation for non-interacting electrons 
with the effective mass m* moving in a trapping potential 
V{r), is given by 



H^.ir) = [-^^ + j = E,'^^{r)■ (1) 

The eigenenergies Ei and eigenfunctions are ob- 

tained using numerical diagonalization. The induced po- 
tential 0ind is then determined from the self-consistent 
integral equation 

0ind(r,c^)= / xo(r',r",w)0tot(r",w)Fc(|r-r'|)dr'd](2) 



Here Vcdr — r'|) is the Coulomb potential and 

Xo{r',r",u) = 

f{E,)-f{E,) 



E 



Ei — Ej — hid ^ i"f 



**(r')vl/.(r")**(r")*,(r') (3) 



is the non-local density-density response function, f{Ei) 
is the Fermi distribution function and 7 is the level broad- 
ening constant. Here 0tot(r, 1^) = </>oxt(r,w) +0ind(r,w) 
is the self-consistent total potential. The external field is 
assumed to be harmonic with frequency w, linearly po- 
larized, and with the wavelength much larger than the 
characteristic system's size, therefore E^xt does not de- 
pend on r, (f>ind(j) = — |Eoxt|a;. The integral equation 
Eq.(2) was discretized on a real-space cubic mesh with 
the lattice constant A, and the resulting system of linear 
equations was solved numerically. In our simulations we 
used the discretization length A = 1.8nm and we set the 
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FIG. 1. The ground state electron density for a nanostructure 
with the trapping potential having the geometry of the Hilbert 
curve. The size of the system L = 60 nm. The number of electrons 
N = 50. 
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FIG. 2. The ground state electron density for a nanostructure with 
the trapping potential having the geometry of the Hilbert curve. 
The size of the system L = 240 nm. The number of electrons 
N = 200. 



effective electron mass to m* — 0.067mo that corresponds 
to doped to the concentration lO^^cm"'^ GaAs TS]. The 
natural energy scale Eq is defined by Eq = fi,^/(2m*A^) 
and Eq « 174 meV in this paper. The level broaden- 
ing constant is set to 7 = 1 x 10~^ Eq. In our simu- 
lations we assumed zero temperature T — OK. The di- 
electric response is quantified using the total energy of 
the local field inside and outside of the nanostructure 
^dr. normalized by the energy of the 
dr. 



W = Jy |i?tot(r) 
applied field Wq = Jy |SextP 

First we calculated the ground state electron density 
for the trapping potential having a constant depth d and 
a geometry of the Hilbert curve. We assumed 50 electrons 
in the system. The potential depth is set to d = 5Eo and 
has geometry the Hilbert curve (Fig.l), The size of the 
system 60nmx60nmx7.2nm. We have also considered a 
larger system with the same electron density. In Fig. 2 we 
plot the ground state electron density assuming N = 200 
electrons in 240nmx240nmx7.2nm. system (please see 
Fig 2). The potential depth is again set to d = 5Eo. 

We calculated the dielectric response of the fractal sys- 
tem shown in Fig. 1 for two electron densities, assuming 
= 25 and A^ = 50 electrons in the system. In both 
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FIG. 3. The normalised local field energy W/Wq near the nanos- 
tructure shown in Fig.l as a function of the frequency of the applied 
field uj for different number of electrons in the system: N = 25 
(solid line), A'^ = 50 (dashed line). 
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FIG. 4. The induced charge density (in arb. units) at the resonant 
frequency cui. The direction of the external field Ecxt('.j) points 
along the y axis. The size of the system L = 60 nm. The dots 
mark the shape of the trapping potential. 



cases the response has a complex dependence on the fre- 
quency of the external field that is consistent with the 
complex geometry of the nanostructure. In the case of 
TV = 25 the normalized energy of the local field W/Wq 
in the system has a distinct maximum at almost zero 
frequency (Fig. 3, solid line). In the case of = 50 the 
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FIG. 5. The induced charge density (in arb. units) at the resonant 
frequency ui2- The direction of the external field Eext('.j) points 
along the y axis. The size of the system L = 60 nm. The dots 
mark the shape of the trapping potential. 
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FIG. 6. The local field intensity logj^g(|£'tot(i")P) at the reso- 
nant frequency u)\, corresponding to the induced charge distribu- 
tion shown in Fig. 4. The direction of the external field Ecxt('.j) 
points along the y axis. The size of the system L = 60 nm. The 
dots mark the shape of the trapping potential. 
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FIG. 7. The local field intensity logj^g(|Etot(i')P) at the reso- 
nant frequency UJ2, corresponding to the induced charge distribu- 
tion shown in Fig. 5. The direction of the external field Eext('i^) 
points along the y axis. The size of the system L = 60 nm. The 
dots mark the shape of the trapping potential. 



response of the system has two clear maxima as a func- 
tion of the frequency of the applied field, one peak is at 
hwi « 0.04ii^o, and the second is at fiw2 « Q.1%Eq (Fig. 3, 
dashed line). We have also calculated the dielectric re- 
sponse for higher electron densities: N = 100, 250 (not 
shown). Higher carrier concentration leads to a shorter 
electron screening length and reduced quantum delocal- 
ization effects [12|. As a result, the dielectric response 
in systems with higher electron densities is more close to 
macroscopic (bulk) systems. The local field enhancement 
is localized near "corners" of the Hilbert curve. The re- 
sponse has many resonances, which are extended over a 
large frequency scale, which is consistent with the com- 
plex self-similar geometry of the nanostructure. We have 
investigated the induced charge spatial distribution at 
resonance frequencies for N = 50. At uji the induced 
charge density has a two dimensional profile: the nanos- 
tructure is divided into upper and lower halves having the 
induced charge densities of the opposite signs (Fig. 4). At 
the frequency of the applied field a;2 the induced charge 
spatial distribution has a pronounced one dimensional 
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FIG. 8. The normalized local field energy W/Wq as a function 
of the frequency of the external field u} for different depth of the 
trapping potential d. Dotted line corresponds to d = 0, solid thin 
line: d = IEq, dashed-dotted line d = 2Eq, dashed line d = 4i?0i 
solid thick line d = SEq- 



electrons are trapped in a quasi-one dimensional poten- 
tial with the Hilbert curve geometry, one can observe 
either two dimensional or one dimensional excitations 
in the system, depending on the excitation frequency. 
We have also studied how the plasmon response depends 
on the depth of the trapping potential, which effectively 
controls the smooth transition of the system geometry 
from quasi-two dimensional to quasi-one dimensional. 
We have found that while the average electron density 
increases with the increasing of the depth of the trap- 
ping potential, the geometry of the Hilbert curve results 
in the overall red shift in the plasmon frequency. This 
opens a broad possibility in control of the resonance prop- 
erties of such systems, making them a robust and repro- 
ducible Raman active nano-structured surfaces for the 
SERS, and other nanoplasmonic applications. 



distribution along the path of the Hilbert curve (Fig. 5). 
The local field intensities produced by these two charge 
distributions are plotted in Figs. (6,7). One can see that 
in the case of quasi one dimensional excitation (Fig. 7) 
the spatial distribution of " hot spots " is much more uni- 
form and has less inactive area, that beneficial for more 
reliable SERS measurements. Note that in Figs. (6,7) we 
used a logarithmic scale. 

We have also calculated the dielectric response of the 
larger nanostructure shown in Fig. 2. for different depth 
of the trapping potential d = 0,1, 2, 4, 5Eq. Based on the 
results of simulations shown in Fig. 3, the maximum of 
the response of a larger fractal-like nanostructure should 
be near w = 0. In the case of d = (finite thin film) the 
system has its response maximum near the surface plas- 
mon resonance (ujpi = uihuik/V^), where bulk plasmon 
frequency luI^^ 



Ik 



— 



:, here p is the carrier concentra- 
tion, e is the electron charge, eg is the vacuum permittiv- 
ity, and e w 11.1 is the high-frequency dielectric constant 
of GaAs [18| . With the increasing of the depth d of the 



trapping potential one naively expects that the resonant 
frequency will shift to a higher frequency, since the same 
amount of electrons will effectively occupy a smaller vol- 
ume, that will correspond to a higher electron density and 
higher resonant frequency. However, the resonance fre- 
quency of the Hilbert geometry has moved to lower values 
(please, see Fig. 8), and at d = 5Eo the resonance becomes 
very close to w — 0. Note that the overall magnitude of 
the normalized local field is decreased with the increase 
of the depth d. One may attribute this to a reduction of 
the screening in quasi-one dimensional systems. As a re- 
sult, the amplitude of the collective (plasmon) resonance 
in the quasi one dimensional system should be less than 
in the quasi-two dimensional system. 

In summary, we have studied the dielectric properties 
of nanoscale systems with the geometry of the Hilbert 
space filling curve. In particular, we found that while 
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